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New Insights on the QSO Radio-Loud/Radio-Quiet Dichotomy: 
SDSS Spectra in the Context of the 4D Eigenvector! Parameter 
Space 
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ABSTRACT 

We search for a dichotomy /bimodality between Radio Loud (RL) and Radio Quiet (RQ) 
Type 1 Active Galactic Nuclei (AGN). We examine several samples of SDSS QSOs with high 
S/N optical spectra and matching FIRST/NVSS radio observations. We use the radio data 
to identify the weakest RL sources with FRII structure to define a RL/RQ boundary which 
corresponds to log Li 4(3ff2=3L6 ergs s^^ Hz^^. We measure properties of broad line H/3 
and Fell emission to define the optical plane of a 4DE1 spectroscopic diagnostic space. The 
RL quasars occupy a much more restricted domain in this optical plane compared to the RQ 
sources, which a 2D Kolmogorov-Smirnov test finds to be highly significant. This tells us that 
the range of BLR kinematics and structure for RL sources is more restricted than for the RQ 
QSOs, which supports the notion of dichotomy. FRII and CD RL sources also show significant 
4DE1 domain differences that likely reflect differences in line of sight orientation (inclined 
vs. face-on respectively) for these two classes. The possibility of a distinct Radio Intermediate 
(RI) population between RQ and RL source is disfavored because a 4DE1 diagnostic space 
comparison shows no difference between RI and RQ sources. We show that searches for 
dichotomy in radio vs. bolometric luminosity diagrams will yield ambiguous results mainly 
because in a reasonably complete sample the radio brightest RQ sources will be numerous 
enough to blur the gap between RQ and RL sources. Within resolution constraints of NVSS 
and FIRST we find no FRI sources among the broad line quasar population. 

Key words: galaxies; active, (galaxies;) quasars; emission fines, (galaxies;) quasars; general 



1 INTRODUCTION 

A much debated problem in AGN studies involves the possibility of 
a real physical dichotomy between radio-loud (RL) and radio-quiet 
(RQ) QSOs. The low fraction o f RL sources - on average ^ 5-25% 
(e.g. Kellermann et al. 1989; Padovani 1993!; iKellermarm et al.l 



1994 ; Jiang et al. 2007) depending on the adopted definition of 
radio-loud ness - and its dependence on redshift and optical lumi- 
nosity (e.g.'Peacock et al.'l98^;'Miller et al.'l990';'Visnovskv et al.' 
;992; Padovani 1993; Hooper et al. 1995; Goldschmidt et al. 1999; 
Jiang et al . 2007) add to the difficulty of defining statistically mean- 
ingful samples of QSOs with which to identify potentially bimodal 
properties. 

Another complication is introduced by the fact that some 
good fraction of RQ sources share common properties with the 
RL quasars; for example: a) about 30-40% of RQ QSOs are spec- 
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troscopically similar to RL (e.g. ISulentic et al.ll2000al and present 
study) and b) both QSO types are capable of producing radio 
jets. RQ jets typic ally extend over scales of a few parsecs up 
to kiloparsecs (e.g. iBlundell & Beaslv 1998; Kukula etal .111998; 
lUlvestad et all |2005| ; iLeipski et aU .2006.fl and RL much larger 
scales with highe r radio power (e.g. iRawlings & Saunderslll99ll ; 
iMiller et al.lll993l) . Potential bimodal properties might be hiding 
behind such similarities. 

The very definition of radio-loudness is rather "loose" with 
continued disagreement over the empirical RL/RQ boundary. Over 
the last few decades a couple of possible boundary criteria have 
been proposed based on: Criterion 1 - radio power JMiller et al.l 
[l99G) and Criterion 2 - radio/optical flux density ratio. Criterion 2 
involves the much used Kellermann factor Rjc (radio flux density at 



^ However lUIvestad et alj 120051) report deep lower frequency VLB I ob- 
servations of several RQ objects studied by Blundell & B easlvl U998h and 
do not confirm the presence of jet-related structure. 



2 Zamfir et al. 



6cm normalized to B-band flux density'). iKellermann et alj jl989l) 
suggested Kk = 10 for the RL-RQ boundary in the Palomar-Green 
(PG) sam ple of QSOs (Schmid t & Green 111 9831 : Icireen et alJI 19861 : 
iBoroson & Green 1992) and many studies have adopted this value. 
Others have suggested that different nominal Rk limits for radio 
steep- and flat-spectrum sources v^ould eliminate the confusion in- 
troduc ed by a fixed value of 10 JFalcke et alJI 199631) . ISikora et al.l 
( 12007 ') propose another kind of quantitative distinction for RL 
AGNs based on Rk as a function of Eddington ratio (see their 
section 4.1). This definition might be relevant if one includes FRI 
sources and LINERs that do not show broad lines. 

There are several different surrogate definitions of Rk in 
literature involving radio measures at various frequencies and 
optical (B-band or i-band), UV or even X - ray measures (e.g . 
Kellermann etahl 1 19891: IStocke et alj Il992l: llvezie et alj 1 2002 : 



Terashima & Wilsonl l2003l : Ijester et al.l 12005 J : IWang et all boOC ; 
Jiang et al. 2007|), which obviously complicates the comparison of 
different studies. 

It is still unclear whether one of the two crite ria is more phys- 
ically significant. Studies like iMiller et alj d 19901) promote the ra- 
dio power as a more fundamental discriminator, while others argue 
in favor of the second criterion, which relates the radio properties 
to other regimes of energy output. Moreover, while some galactic 
nuclei qualify as radio-loud based on one boundary c riterion they 
fail to do so when using the other. A good example dHo & Peng I 
I2OOII) involves a sample of bright Seyfert nuclei where as many 
as ~60% of the sources are RL using Rk > 10 (adopting a "nu- 
clear" radio-optical ratio), but only o ne would satisfy th e condition 
of L(6cm) > 10^^ W Hz~^ s r~^ dMiller et al.lll99(]h . The con- 
clusions in IHo & Peng I ( I2OOII) are provocative in terms of both 
the radio loud fraction and the fact that many of their so-called 
"radio-loud" would be hosted by spiral galaxies, a rather differ- 
ent resul t compared to more l uminous samples anal y zed in studies 
like e.g. [Taylor et alj ( Il996l) , iMcLure et all ( Il999l) , lOunlop et af 




ilOOJ) explains (based on the results of 



Xu et al 



that Rk is Luminosity-dependent and one should rather use 



Rxoc L~ to separate RLs from RQs. The Rk=10 suggested by 
IKellermann et al.l dl989l) for luminous samples (Ms ^-26) is not 
a val id choice for RL boundary for low luminosity samples (e.g. 
IHo & Peng 2001). Nonetheless, Rk retains its heuristic value be- 
cause it offers a scaling relation between nonthermal and thermal 
mechanisms at work in AGN. After all, for theoretical models of 
accretion disk it is preferable to use dimensionless quantities like 
the Eddington-scaled luminosity and accretion rate (although there 
may be different scaling relations for jets, d isk, corona luminosit y 
with accretion rate and black hole mass, e.g. lKording et al.l2006ah . 

Thus, we face the problem of labeling objects differently de- 
pending on the adopted definition of boundary criterion, which 
complicates the integration of different results into a more general 
picture. A further problem involves the combination of radio and 
optical flux measures (the latter being susceptible to internal ex- 
tinction) that can introduce serious selection effects and biases at 
different redshifts thus making Rk a problematic radio-loudness 
indicator for statistical purposes. These aspects reenforce the ne- 
cessity of alternative approaches toward a consistent definition of 
radio-loud. We need a more standardized definition. 

Different studies over the last decade report contradictory re- 
sults regarding the question of a bimodal distribution of QSOs 
in terms of radio-loudness. Rece nt SDSS-based studies (e.g. 
Ilvezic et alj2002l : IWhite et al.l2007h defend the reality of bimodal- 
ity for QSO distribution using histograms of radio/optical(UV) ra- 
tios. The latter study shows a significant dip at Rx~30-40 (see 



their Figure 15). They employ image stacking to lower the de- 
tection limit of FIRST to nano-Jy levels and their final sample 
includes over 41000 sources. A bimodal distribution is found by 
iLiu et alj 12006) in terms of Rk corrected for orientation (although 
with a heterogeneous sample). Their Figure 9 shows a cutoff in 
the RL population in the range logRx = 1.5-2.0. There are at 
the same time many studies that question the reality of bimodal 
ity (JFa lcke et al.' 1996a; White et al. 2000; Brotherton et al. 200 J 
Lacvetal. 2001; Cirasualo et al. 2003a ; Cirasoulo et al. 2003 



Falcke et al 



1 1996iJ) . iLacv et alj ( l200ll) . 



Brotherton et al 



lai- 
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i 200lh 



propose a population of Radio Intermediate sources (RI) that might 
bridge the gap between RL and RQ. Clearly, the RL/RQ problem 
is far from resolved reflected in the lack of consensus on how to 
consistently define a radio-loud sample or prove the existence of a 
physical dichotomy. 

Even if the distribution of radio-loudness measures for a sam- 
ple of RL and RQ sources may not exhibit bimodality, we showed 
that they may represen t two distinct class es of AGN, based on 
spectroscopic measures dSulentic et alj20030 . The present study at- 
tempts to provide more robust empirical support to this alternative 
approach. This paper also shows that the picture of "dichotomy" 
is significantly distorted by mixing bright and faint QSO samples 
when we are flux-limited in b oth optical and radio regimes. 

We recently considered dSulentic et"^l2003l) a third RL/RQ 
boundary criterion (Criterion 3) based on the classical radio mor- 
phology. Double-lobe FRII morphology ( Fanaroff & Riley 1 1 974h 
is the most common type observed in broad line emitting RL 
quasars. FRI morpholog y is very rare among broad -line AGN, 
e.g. 3C120, E1821+64 3 dBlund ell & Raw lingj|200lh and SDSS 
J104022.79+444936.7 dHevwood et al. 200^!The latter reference 
suggests that FRI morphology may become more common among 
broad-line quas ars beyond Z'-^l.O. In a simple orientation unifi- 
cation scenario dUrrv & Padovanill 19951) core-dominated (CD) RL 
sources are interpreted as FRII sources viewed with radio jet axis 
aligned close to our line of sight. If this scenario is valid then the 
CD counterparts of any FRII population will be on a verage more ra- 



dio luminous due to re l ativistic boosting effects (e.g. Orr & Brownd 
1 19821 : IScheueJ 1 19871 ; iBarthell 1 19891 : 1 Jackson & Walll 1 19991) . We 



adopted the weakest (assumed unboosted) FRII sources in our sam- 
ple to define the lower boundary of the RL phenomenon. This 
also allowed us to redefine the boundary in terms of Criterion 1 
(logP6cm~32.0 erg s"^ Hz~^) and Criterion 2 {Rk ~ 70). 

In the present study we reiterate the idea that a robust def- 
inition of a radio-loud quasar can be formulated only if radio- 
morphology is taken into account. More specifically the FRIIs 
should be considered the parent population of RL quasars (e.g. 
Orr & Browne! q^. ischeueij 1 19871 : IBarthell 1 19891 : (Tavlor et alj 
_1996; Jackson & Wall 1999(). This is also supported by the recent 
confirmation that the shape of the observed luminosity function of 
FRII radio galaxies is the same as the intrinsic luminosity func- 
tion of RL quasars dLiu & Zhang||2007l) . In other words, the ra- 
dio weakest FRII structures should dictate the empirical bound- 
ary between RL quasars and the rest of the QSO population. 
We therefore defined a sample of RL quasars using the radio lu- 
minosity coupled with the radio morphology in order to avoid 
the perviously discussed prob lems associated with R^ (see also 
IWadadekar & Kembhavilll999l) . 

We have been ex plori ng a 4D para me ter sp ac e (4DE 1 ; 
iMarziani et all l200ll. l2003j|bl : 
serves as a spectroscopic uni- 
fier/discriminator for all broad emission line AGNs (Type 1). 
Our principal parameters involve measures of: 1) full width at half 
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iSulentic et al 



20()0al lbl: 
20071) that 



SDSS spectra in 4DE1: Radio-Loud/Radio-Quiet Dichotomy 3 



maximum of broad H/9 (FWHM H/9), 2) equivalent widtii ratio 
of optical Fell (A4570A blend) and broad H/3, RFe//=W(FeII 
A4570A)AV(H/3), 3) the soft X-ray photon index (Tsoft) and 4) 
CIVAI549A broad line profile velo city displacement at h alf maxi- 
mum, c(l/2). The "Introduction" of ISulentic et alj j2007h explains 
how this parameter space evolved from various pioneering works. 
One enormous advantage of this parameter space for mulation is its 
weak or abse nt dependence on source luminosity l iSulentic et alJ 
UOOOal . 120040 . Armed with our improved definition of the lower 
boundary for RL activity in quasars we compared their 4DE1 
properties with RQ sources. We found that most RL sources show 
a much restricted domain space occupation within the optical 
parameter plane of 4DE1 (FWHM H/3 vs. Rpeii) compared to the 
RQ majority tSulentic et aL.2003.) . 

Although the results presented in 2003 provided compelling 
support for RL-RQ bimodality, the adopted sample was rather het- 
erogeneous, incomplete and included many sources with measures 
derived from marginal S/N spectra. With the advent of the SDSS 
database it becomes possible to select a large and much more com- 
plete (~ 90%]J sample of AGN with uniformly high resolution 
and S/N optical spectra. A further advantage involves the larger 
wavelength interval 3800-9200A sampled by SDSS. This is com- 
plemented by uniform radio survey data from FIRSTlj designed to 
match the SDSS sky coverage and NVSS survejQ 

The value of studying the radio-loud phenomenon within the 
4DE1 context is at least twofold: 1) it compares RL and RQ sources 
in a parameter space defined b y measures with no ob vious depen- 
dence on the radio properties JMarziani et al.ll2003bl) and 2) it al- 
lows us to make predictions about the probability of radio loudness 
for any population of QSOs with specific optical (or UV) spectro- 
scopic properties. 

The paper is organized as follows: § 2 presents the sample se- 
lection and the RL definition based on radio morphology and radio 
luminosity. § 3 presents the 4DE1 optical measures. § 4 includes 
differences in RQ, RI and RL source occupation within 4DE1. § 5 
offers a discussion on the RL/RQ dichotomy based on L-dependent 
diagrams. In § 6 and § 7 we discuss the fraction of RL quasars and 
the probability of radio-loudness within the 4DE1 optical plane. 
The last two sections are dedicated to discussions and conclusions. 
Throughout this paper we use Ho = 70 km s^^ Mpc^^, Qm = 0.3 
and ^A = 0.7. 



2 DEFINING A POPULATION OF RADIO-LOUD 
QUASARS 

We consider any AGN that shows broad (Balmer) emission lines 
as "QSO" regardless of its intrinsic luminosity (or absolute magni- 
tude). This is why we generated our own sample of SDSS QSOs 
(from Data Release 5 of SDSS: lAdelman-McCarthv et alj|2007l) 
instead of extracting it from the vetted catalog of ISchneider et al.l 
( 120071) that is limited to sources with absolute magnitude Mi < - 
22.0. Sample size is driven by the following goals: i) a sample of 



■^ As defined within the SDSS project the completeness was estimated in 
two ways: by checking how many of the previously known QSOs are re- 
co vered and evaluating th e ou tput of target select i on for simulated quasars, 
see lRichards et alj )2002h and lVanden Berk et alj J2005l) 
^ Faint Images of the Radi o Sky at Twenty- centimeters (FIRST)- 
|http://sundog. stsci.edu/ see also Becker et al. 1 1993) 
"^ NRAO V LA Sky Survey (NVSS) - http://www. cv.nrao.edu/nvss/ 1 see 
also lCondon et al.l < 19981) 



high quality spectra suitable for 4DE1 spectroscopic analysis, ii) 
as complete as possible sample of RL quasars, iii) a large enough 
sample of RQ quasars with reliable 4DE1 measures so that we can 
define the RQ zone of occupation, iv) a representative sample of 
so-called RI sources and v) source-by-source evaluation to avoid 
the pitfalls (e.g. radio/optical misidentifications, misclassifications) 
of automated processing. Our approach is based on careful exam- 
ination of each optical spectrum in order to confirm the presence 
of broad lines. All FIRST/NVSS (20cm/1.4GHz) radio maps were 
visually examined in order to evaluate radio morphology, resolve 
ambiguous cases and obtain the correct integrated (total) radio flux 
density. 

OPTICAL SELECTION: We restricted source selection to 
z < 0.7 so that H/3 and adjacent spectral regions of interest 
(used to define the underlying continuum) would be accessi- 
ble. We selected our sample in several steps. Step 1 selected all 
SDSS DR5 quasars with psf g< 17.0 (n=333 QSOs with 34 RL). 
This selection was motivated by the need for high S/N spec- 
tra from which 4DE1 parameters could be reliably measured. 
Step 2 extended this limit to psf g=17.5 (n=806 QSOs with 76 
RL). This extension was motivated by the desire to increase the 
RL sample. Our fir st two steps are based on a rather blue fil- 
ter, close to BQS JSchmidt & Greeii] Il983l: loreenet alj [198^ : 



BOS 



lBoroson& Greed Il992h . Although Ijester etal.l ( l2005al) find no 
radio-related incompleteness for BQS-like selected samples, they 
point out that the apparently large fraction of RL in BQS survey is 
related to its rather blue filter (B-band) selection (see § 6). Having 
this in mind, we considered also a step 3 aimed to define a RL sam- 
ple considering all QSOs brighter than psf i=17.5 (n=1656 QSOs 
with 91 RL). In Table 1 we explain which objects have been con- 
sidered for the spectroscopic analysis. 

We offer in Figure 1 the sources redshift distributions result- 
ing from the two selections (based on psf g < 17.5 and based on psf 
i < 17.5). Selection using the g-band magnitude limit yields a sam- 
ple with much more uniform redshift distribution. Not surprisingly 
i-band selection (much more complete) favors redder local QSOs 
(50% have z < 0.15; 60% have z < 0.20 and 70% have z< 0.25), 
more strongly affected by host galaxy contamination and/or by the 
presence of Ha line inside i-band. There is obviously a large over- 
lap between the samples. The census presented here also includes 
the objects that show broad lines narrower than 1000 km s^^ and 
are labeled "Galaxy" by the SDSS spectroscopic pipeline (see the 
appendix and see also Table 2). Such sources can be safely included 
under the generic umbrella of "QSO". However, these objects were 
not included in our spectroscopically-processed sample for the fol- 
lowing reasons: 1) limitation of our tem plate used to extract the F ell 
lines from our spectra (same as used in lBoroson & Greenlll992l) or 
2) very red continua, extreme Balmer decrement, significant galaxy 
contamination and dramatically different Ha and H/3 broad lines. 
We also included all QSO spectra that have been assigned only a 
"Science Primary" index of ([j. 

RADIO SELECTION: We used FIRST combined with NVSS 
to evaluate the integrated radio emission and source structure. For 
FIRST survey the typical rms fluctuations are 0.15 mJy, and the 
resolution is 5". For NVSS survey the rms brightness fluctuations 



5 "The SDSS SCIENCEPRIMARY flag indicates whether a spectrum was 
taken as a noimal science spectrum (SCIENCEPRIMARY = 1) or for an- 
other purpose (SCIENCEPRIMARY= 0). The latter category contains qual- 
ity assurance and calibration spectra, or spectra of objects located outside 
of the nominal survey area." iSchneider et alj2007h 
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Table 1. Samples with z < 0.7 selected for spectral analysis in the Context of 4DE1. 



Apparent Magnitude 



^total ^RL ^radio— detected 



Considered for Spectral Analysis? 



psfg< 17.0 


333 


34 


136 


17.0<psfg< 17.5 


473 


42 


122 


psf i < 17.5 ANDpsfg > 17.5 


850 


19 


206 



all 

all radio detected (either FIRST or NVSS) 

onlyifLi.4GWz > 31.0 erg s^^ Hz^^ (allRIandRL) 



Notes: Every QSO in the samples selected based on g-filter was examined in FIRST/NVSS radio maps. 

The radio properties for the sample of QSOs selected by the third set of criteria (based on i-filter) were determined in several steps: 1) we searched FIRST 
within 15" from the optical position, then 2) all sources without a detection (or not in area covered by FIRST) were searched in NVSS around 80" and finally 
3) all sources that are radio-detected in either survey are examined in detail in order to get the total flux density from all components. 



] 1736 QSO with psf i< 17.5 
1 822 QSO with psf g < 17.5 




Figure 1. The distribution of QSOs apparently brighter than psf g : 
(black) and brighter than psf i = 17.5 (light grey). 
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Figure 2. The distribution of logLi.4(3£f2 radio luminosity - calculated 
based on NVSS integrated flux density - for all F RIl quasars identified in 
our sample plus all FRII sources with z < 0.7 from ldeVries et alj I200g) . 



are 0.45 mjy bearn"^, with a 45 " resolution (see footnotes 3 and 
4 of this paper for the sources of this technical details). Both radio 
maps were compared to avoid missing extended sources that might 
have been attenuated with FIRST. The 45" beam of NVSS yields 
sensitivity to more extended structure and provides radio data for a 
few sources not observed by FIRST. We also wanted to clarify the 
nature of any significant discrepancies between the two surveys for 
specific sources. 

We identified all bona fide FRII structures associated with our 
sample of QSOs. Two sources (SDSSJ075407.96+43I6I0.6 and 
SDSSJ080I3I.97+4736I6.0) show bright cores with detached and 
apparently unrelated satellite sources using FIRST. NVSS shows 
that they are FRII with very large (Mpc-scale) radio FRII struc- 
ture that reveals the satellite sources as associated hotspots. SDSS 
J0I352I.67-004402.0 is an excellent example of false FRII and 
false RL. At the same time not all detected double-lobed sources 
in a sample can be unambiguously classified as FRII. However all 
sources with: 1) low enough redshift , 2) a broad line spectrum 
(Type 1 AGN) and 3) adequate radio resolution show FRII (or hy- 
brid e.g. HYMOR - see lGopal-Krishna & Wiital200(]h morphology. 
We therefore assume that all double lobed RL sources in our sam- 
ple are FRII. All sources with FRII structure are assumed to be RL, 
while the other sources (with core or core-jet radio morphology) are 
considered RL only if they have the radio power Li.4ghz above the 
threshold set by the weakest FRII. 

Comparison of NVSS and FIRST flux densities for each RL 



source reveals: a) NVSS measures are larger than corresponding 
FIRST values for virtually all FRII sources and b) most CD RL 
sources show agreement between FIRST and NVSS measures with 
a scatter of --^±2 mJy, most likely due to variability (see also 
IWang et alj 120060 . In a few extreme cases we see differences of 
up to ^-^iSOO mJy. There is no evidence for a significant popula- 
tion of CD RL sources with excess NVSS flux density that might 
be due to extended structure not seen by FIRST. 

We found n=48 FRIIs brighter than psf g=17.5 OR psf i=17.5. 
NVSS flux density measures were always preferred over FIRST 
(for all sources, not only for FRIIs). The small size of the FRII sam- 
ple motivated us to ad d double-lobe quasar s in our redshift range 
that were identified bv ldeVries et al.l 1 I2OO61) based on DR3. Care- 
ful reexamination of optical spectra and radio maps for this addi- 
tion caused us to eliminate several objects (e.g. no broad emission 
lines, no FRII morphology, no radio detection). The DR3-based 
sample, which added n=67 FRIIs, was not optically constrained, 
and most of the quasars in that subsample showed apparent magni- 
tudes fainter than our g or i-band limits. Their sample is not meant 
as an exhaustive list, but the search algorithm had a reported accu- 
racy of ~98% for identifying FRII sources. 

Figure 2 shows the radio luminosity distribution for the 
48-1-67= 11 5 FRII sources. The weakest bone fide FRII/quasar found 
shows logLi,4GHz = 31.6 (erg s"'^ Hz~^). Thus, this becomes 
our radio-luminosity defined boundary between RL and RQ QSOs 
(Criterion 1). The FRII sources in this sample span three orders of 
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Table 2. Additional samples with z < 0.7 used in this study, but not measured spectroscopically. We require that they all show bona-fide Type 1 QSO spectra. 



Apparent Magnitude 



Category /Type 



Ntota! ^radio-detected N_R£, 



psfg< 17.5 

psfi< 17.5 AND psfg> 17.5 

psf g > 17.5 ANDpsfi > 17.5 

19.0<psfg< 19.5 

19.0 < psf i< 19.5 



"Galaxy"-labeled by SDSS 16 12 

"Galaxy"-labeled by SDSS 81 33 1 

double-lobed (FRII) from deVries et al. (2006) 67 67 67 

"QSO"-labeled by SDSS 4800+ 134 47 

"QSO"-labeled by SDSS 3800+ 80 31 



Notes: For the last two samples listed in the table one can notice the very large T^total ^nd a very small T^radio- 
match within one arcsec. We are not concerne d about completenes s of these subsamples used in § 5 and § 6. 
Our sample already included all sources from ldeVries et alj 420061) that have psf g < 17.5 OR psf i < 17.5. 



-detected- We required a SDSS/FIRST optical 



magnitude in radio luminosity (median logLiAGHz = 32.9 erg s~^ 
Hz~^), although they become relatively rare above logLi,4G//z~ 
33.5 erg s^^ Hz^^. Figure 2 shows a continuous distribution of 
radio powers, with a clear decline in the number of sources to- 
ward our RL/RQ nominal boundary. The shape of the distribution 
suggests that there could be only a very small number of FRIIs 
weaker than our weakest bona-fide FRII quasar. It is beyond the 
scope of the paper to fully examine the true nature of the functional 
form describing the FRII radio-luminosity distribution. However, 
we attempted a comparison with the du al population model for 
the radio-luminosity function proposed bv lWillott et al.l j200lh for 
steep spectrum radio sources. We extrapolated the I.4GHz mea- 
sures to 151 MHz measures (to allow for common grounds with 
that study) two ways: 1) assuming a spectral slope ai, of 0.5 and 
2) usin g the empirical scaling relation suggested bv lKording et al.l 
( 120081) (their equation (5)). Either way, we do not see a decline 
on the lowe r luminosity side as abrupt as presented in Figure 3 of 
Iwillott et al.l ( 120011) . Secondly, the Li5ij\/hz distribution that we 
get peaks about one decade fainter (~ 25.5 WHz~^ sr~^) than their 
model suggests for high luminosity population alone (although 
their low-luminosity population could also include FRII sources). 
While not complete, the sample we explore is large enough to be 
assumed representative of the FRII RL phenomenon in broad line 
emitting quasars withi n z = 0.7. 

A recent study dLu et alj 120071) reported eleven extended 
SDSS radio quasars weaker than logLi.4G//z~31.5 erg s^^ Hz^^ 
(see their Figure 2). The authors kindly provided us the list 
of those sources and we analyzed them one-by-one. Two of 
them are included among our RI and are also listed in Table 3 
(SDSS Jl 10717.77+080438.2 and SDSS JI200I4.08-004638.7). 
One other source is also in our sample as RQ with logLi.4GHz 
= 30.7 erg s"^ Hz"^ (SDSS JI62607.24+335915.2). All other 
sources are unambiguously RQ, in some cases off set a few arcsec- 
onds from the optical quasar and were assumed in lLu et al.l J2007I) 
to be associated with the active nucleus. 

The goal of our work was to search for a dichotomy or gap in 
radio properties between RQ and RL sources. Of co urse, sources 
are found with radio inter mediate properties (e.g. iFalcke et al.l 
Il996al : ISuIentic et alj 120030 . Another goal of this paper was to 
isolate a population of these radio-intermediate (RI) sources and 
try to determine if they form a unique (special) class of quasars. 
Unlike the non-RL/RL boundary we have no clear empirical ba- 
sis for defining a RQ/RI boundary because all RI show CD ra- 
dio morphology. The best that we can do is to isolate a region 
that is most clearly RI. Considering that: a) at our sample red- 
shift limit (z=0.7) the minimum detectable radio luminosity (within 
FIRST/NVSS) is logLiAGHz ^ 31.0 erg s"^ Hz'^ and b) the 
radio-FIR (Far Infrared) correlation spans over five decades in lu- 
minosity and extends up to logLi.4GHz — 31.0 erg s~^ Hz~^ (e.g. 



ICondonll 19921 ; lYun et al.ll200l[) . we see that value as a reasonable 
bound ary between star-formation and AGN-domin ate d radio ac- 
tivity (Sopp & Alexander ' 1 99 ll . iKukula et al J[T9 9^ and 'Ha as et~ai] 
I2OO3I show that RQ QSOs follow the radio-FIR correlation). We 
therefore define RI sources as those with logLi.4GHz in the in- 
terval 31.0-31.6 erg s~^ Hz~^ . This is a much mor e restricted RI 
definition than the one given in lFalcke et al.l jl996ah . 

In brief, a source is considered RL if its hi, 40 Hz radio power 
is at least 31.6 erg s^^ Hz^^. A sources is considered RI if its 
Li.4GHz radio power is at least 31.0 erg s^^ Hz^^, but less than 
31.6 erg s~^ Hz~^. All other sources are RQ. 

Summarizing, we considered the whole sample obtained from 
the combination of psf g < 17.5, psf i < 17.5 with an "OR" logi- 
cal operator. The total number of sources was N=1770 (n=95 RL). 
The sample adopted for spectroscopic evaluation (see Table 1) is 
constructed as follows: 

(i) all RL and RI QSOs that are brighter than either psf g = 17.5 
or psf i = 17.5; 

(ii) all RQ QSOs (radio-detected or undetected) that are brighter 
than psf g = 17.0 and 

(iii) all RQ QSOs that are radio detected and show psf g in the 
range 17.0 - 17.5. 

We visually examined the SDSS spectrum for every 
source and rejected objects without emission lines (e.g. 
SDSS J075445. 67+482350.7), objects without broad lines 
(e.g. SDSS J103900.37+414008.7, SDSS J10445I.72+063548.6) 
or with bad pixels (e.g. SDSS Jl 13109.49+31 1405.5, SDSS 
J145638. 81+442755.2, SDSS J220103. 13-005300.2) that pre- 
vented reliable line measures in the region of interest. We rejected 
one supernova: SDSS Jl 13323.97+550415. 8 - as reported by 
IZhou et alj|2003 . Two FRII quasars (SDSS J0924I4.70+030900.8 
and SDSS J123915.39+53I4I4.6) showed serious host galaxy 
contamination (psf g - psf i > I.O in both cases), the broad 
component of H/9 could not be measured. We also excluded from 
our analysis objects with W(H/9) < 20A, which can be sources in 
a high continuum phase (e.g. SDSS JI50324.77+475829.6) and/or 
very red continua with extreme Balmer decrement, where Ha is 
completely different from Uf3 (e.g. SDSS J004508.65+152542.0). 

We should emphasize that our spectroscopic analysis does not 
include rare objects with extreme KfcII > 2.0 values. Sources 
with extrem ely strong Fell ten d to be very red (u-g > 0.8), strong 
IR emitters ( Lipari et al. I993I ; for a detailed study of such a case 
see lVeron-Cettv e t al. 2006) and are not suitable for the Fell tem- 
plate adopted for this study. We identified three RL sources whose 
spectra show extreme Fell emission (RfcII much larger than 
2.0): SDSS J094927.67+3141I0.0, SDSS JI44733.05+345506.7, 
SDSS JI52350.42+3914 05.2. Our attempt to fit the IZwI-based 
dSoroson & Greenlll992i) template for such objects was unsuccess- 
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ful, thus they are not shown in Figure 4. Such objects require spe- 
cial attention and is beyond the purpose of the present study to 
focus on their nature. The inclusion of such pathological and rel- 
atively rare cases will not affect the conclusions of the present 
study. We were able to make reliable measures in the H/3 region 
for N=477 objects. Our RL/RI sample is complete to 17.5 appar- 
ent magnitude in g- OR i-band. We are confident that the RL+RI 
sample is at least 75% radio-complete because all have a FIRST 
S(1.4GHz)> 1.5 mjy (see Figure 1 of Ijiang et alj|2007|). Our RL 
sample includes n=85 RL quasars (n=46 FRII) and n=59 RI QSOs. 
The remaining n=333 represent our RQ sample which, while in- 
complete, is large enough to be representative for the RQ parent 
population. Our final sample includes sources spanning the extinc- 
tion corrected i-band range -27. 5< Mi<-17.1. We cross-checked 
our entire sample selected based on psf g < 17.5 (806 "QSO" -l- 
16 "Galaxy"- s ee Tables 1 and 2) with the "vetted" QSO catalog 
ISchneider et al . (2007). All sources in our sample with Mi brighter 
than -22.0 are present in the that catalog. All sources in that cat- 
alog satisfying our selection criteria are found in our sample. No 
additional QSO satisfying our selection criteria was found there. 




4000 5000 6000 7000 8000 9000 

Figure 3. An example of a QSO contaminated spectrum. The top panel 
shows the initial SDSS spectrum (black) and a template of a host galaxy 
(red); the bottom panel shows the QSO spectrum after subtracting the spec- 
trum of the host galaxy. 



3 ANALYSIS OF OPTICAL SPECTRA FROM SDSS 

In order to obtain the optical parame ters of the 4DE 1 we fol lowed 
the analysis procedure described in iMarziani et al.l yOOSaHj. An 
underlying power-law fx ~ A" continuum was defined using re- 
gions minimally contaminated by Fell lines, specifically at 4195- 
4215 A and 5700-5800 A; it was decided prior to fitting the Fell 
templa te. We used the IZwl-based template of ISoroson & Greed 
( I1992I) for Fell decontamination. Th is represents an impor tant ad- 
vantage over our own Atlas sample JMarziani et alJuOOSal) where 
the typical wavelength coverage of the spectra was ~1000A, ren- 
dering continuum estimation very uncertain. Our chief goals from 
Fell template fitting are: 1) to obtain W(FeII 4750A) blend and 2) 
to clean up the H/9 region. 

In about 20% of sources we attempted to remove the host 
galaxy contaminatior|j before deriving 4DE1 parameters using 
the library of theo r etical galaxy templates from GALAXE\|j of 
ISruzual & Charlod (120030. A muc h more sophisticated approach 
is proposed by IZhou et al.l J2006r) . All our sources with notice- 
able host galaxy contamination lie within z=0.2. There seems to 
be an apparent contradiction with the results of Vanden Berk et al. 
( 120061) . They use an eigenspectrum decomposition technique and 
report a non-zero host galaxy contamination all the way to z=0.7. 
Our sole purpose was to clean up the spectral region of interest of 
the most prominent absorption lines of a host galaxy. We make no 
attempt to estimate the relative proportion of AGN and host light 
in the spectrum unless clear absorption lines are observable; the 
spectra have good S/N to reveal potential host contaminations. Ev- 
idently, most objects we deal with tend to have bluer continua (see 
Figure 1 and the previous section that explains how we built-up our 
sample) and thus are less affected by host galaxy. It is possible that 



^ We used Image Reduction and Analysis Facility (IRAF), distributed by 
the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under coopera- 
tive agreement with the National Science Foundation - http://iraf.noao.edu/ 
^ When prominent absorptions lines like Mg A5177A and Na A5896A are 
present 
* http://www.cida. ve/'~bruzual/bc20()3 



many of the RQ sources not considered for our spectroscopic reduc- 
tion have a significant host galaxy contribution at redshifts higher 
than 0.18. We are also aware that we severely under-represent the 
population of low z objects that form the "bump" in Figure I. As 
we try to explain in this section, the FWHM H/3 seems to be min- 
imally responsive to this uncertainty in the relative amount of flux 
from a host galaxy. See Figure 3 for an example of the host galaxy 
"removal". In most cases this process turned out to be very effective 
in revealing the QSO emission line spectra. We tested the effect of 
host galaxy subtraction on our spectroscopic measures and found 
that ~ 60% of these sources (that required this step) were seri- 
ously contaminated and the rest moderately/weakly contaminated. 
For the latter, we attempted an estimate of the 4DE1 optical plane 
both before and after the correction. We found that FWHM H/3 
changes randomly within ±5%, but Kpeii was much more sen- 
sitive to this procedure, with a 25-35% systematic change toward 
lower Rfe// values. Our set of contaminated spectra span a wide 
range of 4DE1 parameters values. 

When clear inflections were observed between the broad and 
narrow components of H/3, we did not constrain the width of the 
latter component to be the same as FWHM [OIII]AA 4960, 5008 A. 
No attempt was made to decompose the broad Balmer line into 
components. A spline function was used to fit its global profile. 

We compare the values obtained in the pre sent study for 
FWHM H^ and Rfe// with those measured in IMarziani et al.l 
(2003a) for the n=38 sources in common with that Atlas. We cal- 
culated AFWHM H/3 and ARfe// for every object (of these 38) 
considering in each case the values obtained in the current study 
and those in the Atlas (2003). The mean and median differences 
are a reasonable indicator of the robustness of the 4DE1 parameter 
space in its optical dimensions. We find a scatter of ~ 10-20% for 
FWHM H/3 (no systematic effect) and for Rfe// we report a sys- 
tematic effect of 30-35% toward larger values in our present study. 
This latter effect may be due to a more reliable choice of the con- 
tinuum and/or a higher quality of spectral signal. For Figures 4 and 
5 we conservatively adopt the uncertainties estimated in the Atlas 
(2003), even though the quality of the spectra is clearly better in 
our present sample. We are aware that the Kpeii gap (c.f. Figure 
4) between and O.I is not physical and most likely reflects the 
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difficulty of measuring very low values of W(FeII A4570A). Such 
an effect could be due at least in part to our simple definition the 
optical continuum. We would like to explore more on this issue in 
a future project. 



4 LOCUS OF RL AND RI QUASARS IN THE 4DE1 
OPTICAL PLANE 

The optical plane of 4DE1 provides a powerful diagnostic tool for 
testing whether the RL and non-RL sources are spectroscopically 
different. We earlier proposed the existence of two QSO popula- 
tions A and B separated at F WHM H/3::^ 4000 km s"^ (see the In- 
trod uction and section 3. 2 of ISulentic et afluOOTI : see also section 
3 of lSulentic et alj2000ah . 

The previously reported restricted domain occupation of RL 
sources in 4DE1 space means that we can now ask if: 1) the SDSS 
sample confirms the earlier restricted RL domain occupation and 
2) if RI sources show domain occupation more similar to RL or RQ 
Type 1 AGN. Figure 4 shows the distribution of RL and non-RL 
sources in the optical plane of 4DE1 (non-RL include n=333 RQ -l- 
n=59 RI). Figure 4 clearly confirms a restricted domain space occu- 
pation for RL sources with ~ 78% falling ab ove FWHM H/3=40() 
km s~^ (our so-called population B domain: ISulentic et al.l2000ah : 
91% of FRIIs and 62% of CD RL quasars are in this Population B. 
The horizontal line in Figures 4 and 5 marks this nominal popula- 
tions boundary. RQ QSOs show a much wider domain space oc- 
cupation with more than half (~ 62% of our sample) lying below 
FWHM H/9=4000 km s~^. If 4DEI parameters measure aspects of 
Broad Line Region (BLR) kinematics and geometry then this do- 
main occupation difference is consistent with a physical difference 
between RL and the majority of RQ sources, which supports a RQ- 
RL bimodality. At the very least past work reported that RL sources 
show systematically higher black hole masses and systematically 
lower Eddington ratios than the RQ major i ty (e.g. Marziani et al. 
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l200ll : lBorosonll2002l : lMarziani et alj2003bl : ISulentic et al.ll2006l) 

We performed a 2D Kolmogorov-Smimov (K-S) test in or- 
der to evaluate the significance of: 1) the RQ-RL difference and 2) 
the RL FRIIs - RL CDs difference in domain occupation. Follow- 
ing the routine a vailable in Numerical Recipe j j the K-S procedure 
jPeacockl 1 19831 : iFasano & Franceschinil 119871) divides the optical 
plane into quadrants that maximize the two population difference. 
For test 1) the RQ sample is assumed to represent the parent Type 
1 AGN population and RL the test sample. For test 2) the RL FRII 
sample is assumed to represent the parent population and the RL 
CDs the test sample. The results of the two tests are sunmiarized in 
Table 3. One can notice that the RL/non-RL separation at FWHM 
up = 3875 km s~^ is in reasonable agreement with our previously 
adopted Population A/B boundary at 4000 km s^^ (Figure 4). As 
reported in Table 3, the probability that RL and RQ occupy the 
same spectroscopic domain is very low. Similarly, the second test 
between FRIIs and RL CDs shows that the two samples are very 
distinct in terms of spectroscopic properties. In the former case the 
result is equivalent to saying that it is extremely unlikely that RL 
and non-RL are drawn from the same parent population. In the later 
test, the probability listed in Table 3 could be interpreted in two 
ways: I) the orientation is responsible for the distinct space occu- 
pation for FRIIs and RL CDs or 2) the CD RLs and the FRIIs are 
drawn from distinct populations, in which case the RL CDs (all 
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Figure 4. RL quasars in the optical plane of the 4DE1 pai'ameter space. 
The green arrow indicates the displacement between the median FWHM 
H^ and the median Rfe/z for the FRII and the CD sources. The solid light 
gray symbols are non-RL objects. In the upper right com er are indicated 
the typical la errors, estimated in lMarziani et al.l <2003al) . The red dotted 
lines show the boundaries for the quadrants defined by the 2D K-S test we 
performed as we explained in the text. The black dotted line indicates the 
Population A/B boundary. The vertical axis is truncated at 12000 km s^^ 
for clarity, thus we miss showing five other RL (two CD and three FRIIs 
between 12000-24000 km s"!). 



or most of them) could be interpreted either as beamed RQs or as 
pre-/postcursors of an FRII episode. 

There is also a clear bimodality of RL/RQ in terms of Kpeii- 
(robustly confirmed by the 2D K-S test we performed). In the con- 
text of 4DE1 we focus on FWHM H/3 for several reasons: 1) it is 
a direct measure of the Broad Line Region (BLR) kinematics, 2) 
it can be measured more accurately than Rpeii (see § 3) and 3) 
there are significant Population A/B differences reported over the 
last seven years that are defined in terms of FWHM H/3 alone (see 
Table 5 in Sulentic et al. 2007) . 

Figure 4 then shows a significant displacement between the 
non-RL and RL distributions with most RL lying above FWHM 
H/9=4000 km s~^. We also observe a separation between the 
mean/median position of FRII RL sources and core/core-j et (CD) 
RL so urces (see Table 4), confirming a result from Sulentic et al.l 
( 120031) . This is the first step in estimating the role of source orienta- 
tion in 4DE1. The vector shown in Figure 4 indicates the change 
in median 4DE1 optical parameters between FRII and CD RL 
sources. Orientation-unification scenarios see the latter sources as 
having radio jets aligned to within a few degrees of our line of 
sight. The vector suggests that source orientation strongly influ- 
ences FWHM H/3 measures and, to a lesser degree, Rpeii- Given 
the likelihood that large disk-jet misalignments can occur and that 
radio structure in many RL is highly nonlinear it is surprising how 
large is the FRII-CD median parameter separation. The 10-20% 
of CD and FRII RL with, respectively, very large and very small 
FWHM values could be interpreted as sources where the radio 
structure and the region that emits the broad lines (i.e. accretion 
disk) appear to be misaligned, if one invokes the unified picture for 
AGNs relative to Figure 4. The few RL CDs with very large FWHM 
H/3 may, apparently disconnected from the bulk of RL CDs, may 
be the best candidates for a pre- or post- FRII phase. Sources with 
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Table 3. Two-dimensional Kolmogorov-Smimov tests for RL / non-RL and for RL FRIl / RL CD. 



Samples 

nl-parent and n2-test 


Coordinates of quadrants 
FWHM H/3 (km s^^) ; RfeII 


Probability 
of null hypothesis 


392-non-RL and 85-RL 


3875 ; 0.49 


P~6.2xlO-** 


46-RL FRII and 39-RL CD 


6100; 0.18 


P~9.8xlO-* 
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Figure 5. RI QSO in the optical plane of the 4DE1 parameter space. The 
light gray symbols show the RQ objects and the solid red stars are the 
n=59 RI sources. SDSS J232721.97-I-152437.3 is indicated with an "X" - 
see "Discussion" and also Table 3. The vertical axis is truncated at 12000 
km s^^ for clarity, thus we miss showing two other RI. 



extremely broad Balmer profiles (sometimes double peaked) are 
so rare (we find a handful of such sources with FWHM H/3 in the 
range 12000 - 30000 km s^^) as to defy interpretation as the simple 
tail of a normal QSO FWHM distribution. 

Figure 4 included RI along with the RQ QSOs as non-RL 
QSOs. The tacit assumption was that RI and RQ are the same. Fig- 
ure 5 presents a test of that assumption that is equivalent to that 
performed for RL. Does the previously defined RI sample show 
a distribution in the optical plane of 4DE1 more similar to RL or 
RQ sources? The n=59 RI sources (stars in Figure 5) show no dis- 
tinguishable difference in occupation from the n=333 RQ sample. 
Only 42% of the RI population is found in the Population B do- 
main, comparable to the ~ 37% for the RQ sources. There is there- 
fore no spectroscopic evidence that the RI sources form a special 
class (see Table 4). 



5 CAN WE REVEAL A RL/RQ DICHOTOMY USING 
L-DEPENDENT DIAGRAMS? 

The previous section compared RQ, RL and RI sources in a 
Luminosity-independent context. We now address the problem of a 
RL/RQ dichotomy from a Luminosity-dependent perspective. Fig- 
ure 6 plots source bolometric versus radio luminosity (htoi vs. 
Li.4Gifz)- The radio luminosity is K-corrected ( Hogg 1999) assum- 
ing that fv ~ i^" and a = - 0.5 in the radio regime. The bolomet- 
ric luminosity was estimated from ht,oi — W\Lx, where A = 5400A 
(see a concise discussion on the bolometric correction in section 2.8 
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Figure 6. The distribution of our relatively bright sources in a plane defined 
by the bolometric and radio luminosity. The thi'ee objects marked with an 
X are commented in Table 3 and related discussion in § § 8.1. The radio- 
undetected QSOs are not shown here. 



of lMarziaiii et alj2006l and references therein). The 5400A specific 
luminosity is estimated using the continuum flux in the rest-frame 
spectrum of the QSO. We used dopcor task in IRAF with the the 
appropriate cosmological flux corrections applied when deredshift- 
ing the spectra. Our sample covers about four decades in bolometric 
luminosity with logLtoi ~ 43.0-47.0 erg s~^. 

There are a few important comments about Figure 6: a) we 
find no RL quasar fainter than logLi,oi = 44.3 erg s^^, b) we 
find no CD RL below logL^o! ~ 45.0 erg s~^, c) the RQ sam- 
ple shows a power-law correlation {hboi oc \^iaghz^'^^) (see also 
e.g. lKukulaetalJll99§) or altematiy elv Liaghz oc Lboi°'^^ (the 
same as reported in lWhite et alj2007l considering that hboi oc hopt) 
but d) the RL population however shows a different behavior for 
CDs and FRIIs. The majority of CDs concentrate at high values of 
hboi, while the distribution of FRII quasars shows a rough corre- 
lation parallel to the RQ one. Most of the CDs, especially below 
loghi AGHz — 32.5, make no sense in an orientation unification 
scenario because there is no corresponding weaker FRII popula- 
tion to the left of them from which they could be boosted. As we 
suggested in § 4 based on the 2D K-S test for FRIIs/RL CDs these 
sources could be seen as boosted RQs. 

4DE1 parameters show no obvious dependence on source lu- 
minosity. They also show no dependence on radio luminosity apart 
from a restricted domain space occupation observed only for RL 
sources. If RQ and RL sources belong to the same family then 
we could reasonably expect them to follow the same correlation 
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Table 4. Mean and Median optical spectroscopic measures for the samples used in Figures 4 and 5. 



meanit standard deviation 



median 



FWHMH/3 Rpell FWHMH/3 Rpeii 
(kms-i) (kms-i) 



Relative to Figure 4 


RL (n=85) 

non-RL (n=392 RQ+RI) 


6809±3976 
4016±2569 


0.36±0.29 
0.57±0.38 


5775 
3375 


0.27 
0.50 


RL FRll (n=46) 
RL CD (n=39) 


7673±3733 
5789±4059 


0.26±0.17 
0.48±0.36 


6750 
4418 


0.20 
0.38 


Relative to Figure 5 


RQ (n=333) 
Rl (n=59) 


3852±2114 
4941 ±4230 


0.56±0.35 
0.62±0.52 


3227 
3659 


0.50 
0.50 



between bolometric and radio luminosity. Figure 6 suggests that 
RQ and RL sources show separate correlations. We see a clear 
dichotomy between the populations at our previously determined 
boundary (logLi,4GH3=31.6 erg s~^ Hz~^) as lower limit for RL 
activity. This dichotomy appears to independently confirm our pre- 
vious suggestion, based on 4DE1 occupation, that RQ and RL 
sources are fundamentally different. 

Figure 6 is based upon our bright SDSS sample. Does it in- 
clude all FRII quasars within z=0.7? What is the effect on Figure 6 
of including fainter sources in the same redshift range? Considera- 
tion of these questions can help us understand why such conflicting 
results about a RQ-RL dichotomy /bimodality have been obtained in 
past studies. We have reason to fear that Figure 6 does not tell the 
full story about dichotomy because almost all sources in our sample 
with logL(,o;<45.5 erg s~^ show z<0. 15 while all sources above 
that value show z>0.15 (see Figure 8). We have essentially sam- 
pled the bright end of the low redshift Optical Luminosity Func- 
tion (OLF). On this bright end RL sources are relatively abundant 
and RQ numbers small enough to allow a dichotomy to be seen. 
But we have severely undersampled the faint end of the OLF (e.i 
iBovle et alj[2000l : ICroom et al.ll2004l : iRichards et ai]|2005l. 12006 
In that luminosity range the RQ population is so large that the radio 
bright tail of the RQ distribution might overlap the RL distribution 
effectively quenching any dichotomy. 

The above suggestion can be tested and illustrated by adding 
fainter subsamples of QSOs (see Table 2) to Figure 6, leading to 
Figu re 7, which shows th e following : i) the n=67 sample of FRII 
from ldeVries et alj ( l200q) , fainter than our sample of FRII (see § 2); 
ii) all radio-detected objects that are labeled "Galaxy" by the SDSS 
spectro-pipeline but show Type 1 spectra (no spectroscopic reduc- 
tion was performed on these, as we explained earlier); iii) core ra- 
dio sources (no FRIIs) with psf g in the range 19.0-19.5 and iv) 
core radio sources (no FRIIs) with psf i in the range 19.0-19.5. In 
order to understand some of the the subtle effects that come into 
play at this point Figure 7 should be approached in conjunction 
with Figures 8 and 9. We point out that for Figures 7 and 8 we esti- 
mated the bolometric luminosity following the empirical results of 
iHopkins et al.l l2007h - first estimating the B-band luminosity for 
our sources obtained from psf u and psf g magnitudes (corrected 
for extinction, using the SDSS coefficients). The B-band luminos- 
ity (K-corrected, assuming that f^'^v", where a=-0.5 in the optical 
domain) was obtained from psf u a nd psf g magnitude s using the 
transformation formula proposed bv ljester et alj J2005al) . We make 
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9 RL-FRII -deVriesetal. 2006 

□ CD QSO with psf = 19.0-19.5 

O CD QSO with psf 1 = 19.0-19.5 

• the radio-detected ' ' Galaxy" sample from Table 2 
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Figure 7. The distribution of our relatively bright sources plus several sub- 
samples of QSOs selected as explained in the text. The radio-undetected 
QSOs are not shown here. The vertical dotted lines at 31.0 and 31.2 repre- 
sent the minimum detectable radio luminosity that coiTesponds to a detec- 
tion limit of 0.7 - 1.0 mJy at z=0.7, respectively. 



no attempt to reduce the spectra of these faint objects and therefore 
no 5400A specific luminosity (in the underlying continuum) can be 
estimated and used to esti mate hboi- Th i s is w hy we employ the 
the empirical recipe from iHopkins et al.l OOOTJ) . It is worth men- 
tioning that we examined every optical spectrum and radio image 
in FIRST/NVSS to confirm Type 1 optical and FRII/CD radio sta- 
tus. No attempt was made to estimate a host galaxy contribution for 
these fainter sources. The main results of this study (e.g. Figures 6, 
7 and 8) are not sensitive to the choice of the bolometric correction. 
We also remind the reader that the distribution in Figure 9 reflects 
the heterogeneous construction of the radio sample we investigated, 
as explained in the previous sections. 

In Figure 7 it is now obvious that the optically faint radio core 
sources (the empty squares and diamonds) tend to fill in the region 
logLboi ~ 44.0 - 46.0 and logLi aghz ~ 30.5 - 32.0. The fact that 
such objects are less luminous than logLbo; = 46.0 is not a surprise. 
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Figure 8. The luminosity distribution witli redshift for the various subsam- 
ples of QSOs employed in Figure 7; the symbols are similar to those in Fig- 
ure 7, with the exception that this time all QSOs in our sample (RQ, RL CDs 
and RL FRIIs) are displayed with the same solid lig ht grey symbols. Th e 
crossed-circles denote only the FRII subsample from ldeVries et alj j200a) . 
just as in Figure 7. 
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Figure 9. The radio luminosity distribution with redshift for the various 
subsamples of QSOs employed in Figure 7 and 8; the symbols are identical 
to those in Figure 8. 



They are selected as apparently fainter objects in the same z range 
as the bright sample. It is important noting that these fainter sam- 
ples (iii and iv above) are not represented below z=0. 1 and z=0.2, 
respectively. This is the main cause of their tendency to scatter 
mostly in the region that was previously (in Figure 6) scarcely pop- 
ulated along the radio-luminosity axis. The main conclusion we get 
is that the picture of the dichotomy becomes "blurry " at this point. 
W ith the addition of the fainter FRIIs from IdeVries et al.l 
( 12006 ) we observe the parallel RL FRII sequence more clearly in 
Figure 7. The trend we see for FRII quasars is qualitatively simi- 
lar to that reported by other stud ies for ste ep-spectru m quasars us- 
ing lo wer radio frequencies (e.g. lSerieant e t al. 1998 dWillott et al.l 
119981) or higher radio frequencies jXu et alJ1999l) . The two correla- 



tions likely indica te a different fraction of power output channele d 
into the jets (e.g. iRawlings & SaundersI [19911 : iMiller et al.lll993h . 
The core of RL population covers ~2.5dex in the Lto; space and 
~2dex in the Liaghz space, while the RQ sample distributes over 
^4dex in either measure. The small number of extreme sources 
above logLi,o;~46.5 lie near our high redshift limit and suggest that 
we are seeing the first hints of luminosity evolution in our sample. 



6 ESTIMATING THE RADIO LOUD FRACTION (RLE) 
OE TYPE 1 AGN 

We are now in a position to estimate the fraction Radio Loud 
quasars relative to the total population of QSOs that satisfies our 
redshift and luminosity criteria. Figure 1 shows that an i-band se- 
lection adds a large number of relatively low redshift sources. Bluer 
(g-band) selection provides a more uniform sampling of sources in 
our redshift range. If we consider only sources with psf g < 17.5 
the RLF ~ 9.2% (of n ~ 822 sources, i.e. 806 QSO + 16 "Galaxy" - 
see Tables I and 2). Considering sources selected with psf i < 17.5 
yields a RLF ~ 5.2% (relative to n ~ 1737 objects, i.e. 1656 "QSO" 
-1-81 "Galaxy" - see Tables 1 and 2). If we combine the above con- 
ditions with a logical "OR" operator we get RLF ~ 5.4%, relative 
to n '^ 1770 objects. 

Our results are quite different compared to the RLF estimate 
for the PG sample of QSOs; RLF ~ 17% (appl ying our definition 
of rad io-loudness to the sample of 87 sources in lBoroson & GreenI 
1 1 9921 or directly to the full li st of 96 Palomar-G reen UV-excess 
selected QSOs/Seyferts from ' Green et alj Il986l with z < 0.5). 
Boroson & Green ( 1992) include 50 Population A sources (two of 
them RL) and 37 Population B sources (twelve of them RL). Ob- 
viously, one must also bear in m ind that the absolute ma gnitude 
cut at Ms =-23 in the PG sample JBoroson & Green|[T993) has its 
play into the rather large RLF it includes, considering also our re- 
sults presented in Figure 6; we find no RL in the dimmest decade of 
bolometric luminosity we sampled. This is also consistent with the 
reported dependence of RLF on luminosity (see the "Introduction" 
for r eferences re l ated to this issue). 

IJester et alj f2005a) offered a detailed discussion on possible 
radio-re lated incompleteness in the Palomar-Green Bright Quasar 
Survey JSchmidt & Greenlll983l : ICreen et al.lll986h compared to 
an i-band limited sample from the SDSS. They suggest that the 
rather large RLF in the BQS appears connected to the fact that BQS 
objects, being selected in a B-band flux-limited survey, have rather 
blue continua, objects with bluer continua apparently tend to have 
stronger [OIII] lines, and objects with larger [OIII] lines are more 
likely to be radio-loud. In the present study we also find a RLF 
fraction approximately double when selecting the sample based on 
a bluer filter (g-band) than when selecting bas ed on a redder one 
(i-band). On the other hand IJester et al.l ( l2005ah find no systematic 
radio-related biases by comparing the BQS sample against a BQS- 
like sample selected from SDSS database. More recently, two other 
SDSS/FIRST - based studies (deVries et al. 2006; Lu et al. 200^ 
found that FRII quasars are rather bluer than the radio-compact 
sources, which apparently could be connected to the RL excess in 
BQS survey. Indirectly confirming their conclusion (i.e. corrobo- 
rated with our results that the large majority of FRIIs are members 
of Population B) is the study of Richards et al. (2003), which report 
a systematic narrowing of the H/3 line with increasing redness. 

However, the cause of the RL excess in PG survey is not com- 
pletely clear at this time, considering also the apparently contra- 
dictory conclusions about the optical colors of radio quasars, i.e. 
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RL quasars ha ve been found to be in general redder than RQ QSOs 
(e.g. Brothert onet al.ll200ll : llvezicet aljbood ; IWhite et alJboOTi : 
iLabita et al.ll2o'o7l) ! 



7 4DE1 OPTICAL PLANE AND THE RADIO-LOUDNESS 
PROBABILITY 

4DE1 is a diagnostic tool that could set empirical constraints for 
the theoretical models of AGN physics. Our results show that RL 
quasars prefer a restricted zone of occupation in the optical plane 
of 4DE1 (Figure 4). Our g-band selected sample of n=333 QSOs 
(n=34 RL) is the most complete that we have available for this 
analysis. The fractions of Population A and B sources are ~60% 
and ~40% respectively of which 4-5% and 17% are RL. The situa- 
tion is even more extreme if we correct sources for line-of-sight 
orientation (see Figure 4). In an orientation unification scenario 
many of the Population A RLs (mostly CDs) could simply be face- 
on oriented Population B RLs. FRII sources can be said to define 
the locus of the RL population in the 4DE1 optical plane. If one 
extrapolates the relative proportion of Population A/B sources to 
the whole sample of n=1770 (see the previous section), the radio- 
loudness probability would be 10% and 2% for Population B and 
A respectively. A quasar is approximately 4x more likely to be RL 
if it shows a Population B optical spectrum. The scarcity of RL 
qua sars in the population A domain is not in dispute. Recent stud- 
ies jKomossa et alj|2006r) searching for RL Narrow-Line Seyfert 
1 (NLSyl) sources (extreme population A QSOs with FWHM 
H;3 < 2000 km s~^) find that most so-called RL NLSyl are in 
fact RI QSOs. A very s mall number of bona fide RL NLSyl is 
known at this time (e.g. I zhou & Wand |2002| : IZhou et alj |2003|. 
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8 DISCUSSION 

One of the most fundamental differences among the broad line 
QSOs involves the existence of RQ and RL populations. Do all 
sources pass through a RL phase? Do RL quasars represent in some 
way a physically distinct class? The latter question motivates the 
high level of interest in the possibility of a RQ/RL dichotomy. We 
have shown that plots in terms of radio and bolometric luminosity 
(Figures 6-8) are not the best way to answer it. Samples restricted 
to high Lboj sources show hints of a dichotomy, but more complete 
samples do not. This is due to the numerical imbalance in the two 
populations. A complete sample of QSOs contains so many RQ 
sources that the radio brightest of that population will bridge the 
gap between RQ and RL. 

4DE1 provides a better way to address the problem. A way 
that is independent of the radio and optical luminosity of sources. 
4DE1 suggests that the answer to the latter question may be "yes". 
If 4DE1 parameters measure fundamental aspects of BLR kine- 
matics and geometry then we have evidence that the RL quasars 
(mostly Population B) may be significantly different (§ 4) at a fun- 
damental level from the majority of RQ QSOs which occupy the 
Population A domain. Taken at face value, the 4DE1 optical plane 
suggests necessary (yet not sufficient) empirical constraints for de- 
veloping RL activity. It is important to remember that the RQ-RL 
separation in 4DE1 is not complete. About 60% of RQ QSOs (Pop- 
ulation A) show properties almost never seen in RL sources while 
about 40% of RQ sources are spectroscopically indistinguishable 



from RLs (Population B). In this case Population B RQ sources ap- 
parently have the necessary BLR properties for radio loudness but 
not sufficient to be RL. 

One possibility is that population B RQ (and especially RI) 
sources might be the pre- or post-cursors of the RL phase. How- 
ever, another scenario is that the population B RQ sources oc- 
cupying the RL domain simply reflect the overlap of two un- 
related AGN sequences. The idea of two distinct populations 
suggests that something else is a necessary ingredient in AGN 
physics that manifests/tr iggers radio loudness; i t has been sug- 
gested the BH spin (e.g. 'Wilson & Colbert|[l993 : iModerski et all 
1998; Meier 2001; Volonteri et al. 200l ), the host galaxy m orphol- 
ogy (e.g. Capett i & Balam verde 2006; Sikora et al. 2007) and/or 
its link w ith the nucleus (e.g. .Hamilton et aL.2008.) , the environ- 
ment (e.g lKauffmarm et al.l 120070 . Some more clues could come 
from a n analogy with X-ray binaries (e.g. iMaccarone et alj|2003l ; 
Ijester 2005b; Kording et al. 2006a, b). Those studies suggest that 
RL quasars are in a distinct accretion mode compared to RQ QSOs. 
Moreover, for a better understanding of radio-loudness one should 
also consider the ratio of optical:X-ray emission (i.e. diskxorona 
relative emission)(Kording et alj |2006bl) . Future studies can cer- 
tainly explore these valuable arguments incorporating the empirical 
data presented here. 

The overlap of RL and RQ sources in the Population B domain 
suggests that population A-B distinction may be more physical than 
RQ-RL comparisons (see also Boroson 2002). The two populations 
(A and B) are nominally separated at FWHM H/3=4000 km s" ^ and 
there are many other forms of evid ence that support a boundary 
near this value JSulentic et al.ll2007n . We suggested that this might 
correspond to a critical Eddington ratio (log L/Lfitij ~0.15) where 
the BLR p roperties change rather suddenly (e.g. Sulent ic et al.l 
l2000bl;lMai-ziani et al. 200lll2003bll2006; Sulentic et al. 20071). In 



rece nt study JKellv et al.ll2008l with reference to iBonning et al.l 
1 20070 argue that LboifLEdd ~ 0.3 could indicate some critical 
change in the accretion disk structure. At this time we can only 
add that all objects in our sample that have htoifLEdd larger than 
this value are exclusively part of the Population A while all others 
showing htoif^Edd less than 0.3 are a mixture of Population A and 
B. 

Comparing median FWHM H/3 values for FRII and CD RL 
sources gives 6750 and 4400 km s^^ respectively. This differ- 
ence is interpreted as a manifestation of source orientation. CDs 
viewed as near disk face-on (alternatively jet-aligned) sources show 
FWHM measures that are not dominated by Keplerian motions in 
contrast to FRII quasars. Detection of this FRII-CD difference in 
median FWHM (se e also lRokaki et alj|2003l ; ISulentic et alj|2003l; 
IdeVries et alj 120060 supports BLR models involving a flattened, 
disk-like geometry (for a more det ailed discussion on BLR struc- 
ture and dynamics see section 3.1 in lCoUin et al.l2 006V The results 
on the relative distribution of RQ, as well as FRII and CD RLs, 
in Figure 4 indicate that face-on RL sources contain a significant 
extra (~2-3000 km s~^) component of line-of sight motion that is 
not present in face-on RQ AGN. That may or may not be associated 
with the radio jets. 

8.1 Radio Intermediates 

Armed with evidence that BLR structure in RQ and RL sources 
may be fundamentally different we return to the sources with log 
Li.4GHz~31.0-31.6 erg s~^ Hz~^. They are one, or some combi- 
nation, of the following: a) radio-weakest RLs, b) radio-strongest 
RQs or c) a special class (RI) of sources. We disfavor interpre- 
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tation a) because we find no bona-fide FRII radio morphologies 
among them. We favor option b) because they show CD emission 
(Mke RQs), but cannot be boosted FRIIs. Claims of relativistic jet 
detection (e.g. Blund ell & Beasly 1998) in some of these quasars 
has led to the sug gestion that they might be boosted RQ (option 
b) sou rces (e.g. iMiller et al.lll993l : lFalcke et alj[l9963 . IWang eTal] 
|200g), which would even more effectively bridge the gap between 
the majority of RQ sources and the RL quasars. Since the RI show 
no obvious difference from weaker and unambiguously RQ sources 
we disfavor option c) and again favor option b). 

It has also been proposed that the radio emission in 
RQs (or alternatively, non-RLs) is mostly related to star for- 
mation processes (circumnuclear starbursts, supemovae) (e.g . 
SoDP & AlexandeJll99ll; iTerlevich et alj|l992l : [Miller et al.|[T993 



Colina & Perez-01ealll995[) . This idea appears to be naturally re- 



lated to the fact that RQ sources follow the radio-FIR corre- 
lation, as mentioned earlier. On the other hand, the discovery 
of flat radio spectra, elongated radio cores in non-RL quasars 
and/or high brightness temperatures (e. g. iFalcke et al.l Il996al : 



iBlundell & Beaslvll 19981 : IWang et al.ll2006l : TLeipski et alj|2006l) fa- 
vors the hypothesis that the mechanism of radio emission in non- 
RLs is similar to that of RLs. However, the cl aim of a relativistic 
jet in "RQ" PG1407-I-265 JBlundell et ai]|2003l) involves an unam- 
biguously RL quasar by our definition (log Li.4G jfz=32.5 erg s~^ 
Hz~^). Another deep search JUlvestad et al.l2005n failed to confirm 
some of the other detections in lBlundell & Beaslvl (1998) indicat- 
ing that the frequency of occurrence of weak jets in RQ quasars is 
still uncertain. Some of the claims involve AGN that do not show 
Type 1 spectra (e.g iFalcke et al.l |2000| : iNagar et all I2OO0I . I2OOII : 
lOarcia-Baretto et alj|2002r) . which are not considered here. 

All RI sources show core (or core-jet) morphology leading us 
to conclude that their radio emission may be fundamentally differ- 
ent from the classical RL sources (option c?). They could be frus- 
trated jets that face too much resistance from the ambient medium, 
thus failing to manifest as large scale FRII structures. They cannot 
be boosted classical RLs unless we have missed a significant un- 
boosted FRII population which would presumably lie in the zone 
of our RI sample. Our comparison of NVSS and FIRST fluxes for 
the RI sample does not allow us to rule out the existence of hid- 
den FRIIs. A total of n=37 RI were detected in both radio surveys 
and we find about 18 QSOs with an NVSS flux density excess in 
the 10-50% range. Some of these involves unrelated point sources 
within the NVSS beam. A complementary approach is to look for 
unusual radio structures in the FIRST and NVSS maps. 

Table 3 summarizes the properties of quasars showing unusual 
radio structure. Three of these sources likely involve weak FRII 
structure while the rest show no hint of it. These three sources are 
marked with an X in Figure 6. Their location along the trend de- 
scribed by the bona-fide FRIIs increases the likelihood that they 
may be very weak FRIIs. 

Since they are not classical RL sources some of the RI might 
be pre- and/or post-cursors of the classical RL phenomenon. Sup- 
port for this interpretation might come from observations showing 
a flat or curved radio spectra JFalcke et al.lll996al : lb'Dealll998h . 
Sources in this restricted regime merit multifrequency radio mea- 
sures. The most interesting RI source in this context involves 
SDSS 1232721.97-^152437.3 (Table 3) that is indicated with an 
"X" in Figure 5. It is RI based on an integrated radio luminosity 
logLi.4G_ff2=31.2 erg s^^ Hz~^. The NVSS radio map shows very 
extended weak lobes and a strong core (NVSS core/lobe'~3.4). It 
was not observed by FIRST but it is unlikely that FIRST would 
have detected the very extended lobes. The low surface brightness 



in these lobes may be the signature of a past episode of RL ac- 
tivity with the bright core possibly indicating a renewed phase of 
radio activity. We may be observing this source between radio out- 
bursts when old decaying lobes can still be detected. A possibly 
related RL analog involves SDSS 1110538.99-^020257.4. FIRST 
detects a strong core elongated ~45° to the direction of very low 
surface brightness lobes, only one bright enough to be listed in the 
FIRST source catalog. This is likely another example of a two- 
phase RL with very old and very young lobe structures. Without 
doubt a few RI and/or RQ involve weak lobes but there is no evi- 
dence for a large population that could boost many CD sources in 
the range logLi.4G_frz=3I.0-31.6 erg s~^ Hz~^. Or one could ex- 
tend this range to include all of the so-called RL CD s fainter than 
logLi 4Ggz~32.5. Deeper m aps (e.g. PG 1 309-1-355 in lFalcke et al.l 
Il996bl : luivestad et al.l2005l) have failed to turn up weak lobe struc- 
tures. 

The RL boundary is also s imilar to the FRI/FR II break at ~ 
logLi.4Girz=32 erg s^^ Hz^^ ( Owen & Laindl 19891) . It is impor- 
tant to point out that the FRI/FRII break is not a sh arp one and may 
be a function of optical as we ll as radio luminosity JOwen & White! 
[1991; Ledlow & Owen 199o). The recent deep radio survey of RQ 
quasars (Leipski et al. 2006) revealed several RQs with elongated, 
complex and even double sided structures. One source (out of 14) 
shows (VLA B-array) structure (PG0026-I-I29) reminiscent of FRI 
morphology but on a very small (~1.5 kpc) scale and is 1.5 orders 
of magnitude less radio luminous than the weakest FRII in our RL 
sample leading us to conclude that this level of activity is unrelated 
to the classical RL activity. FRI sources are essentially absent from 
our type 1 QSO sample. We found no FRI structures on FIRST and 
NVSS maps for any of our RL/RI/RQ sources (Table 3). 



8.2 Biases? 

Is there a chance that we missed some RI, or especially RL QSOs, 
in Population A? These could expand the RL domain in Figure 4 
lessening the strong RQ-RL difference that we found. As pointed 
out before we should have detected all RI and RL sources in our op- 
tically selected SDSS samples. We consider four possible sources 
of bias in this study. 

a) It is known that SDSS is biased against very narrow 
broad emission line QSOs (often called Narrow Line Seyfert 
rs=NLSyIs) because at least one line with FWHM>IOOO kms^^ 
is required to be assigned QSO type. We attempted to identify 
all such extreme NLSyl sources within our redshift and magni- 
tude limits that are assigned "Galaxy" type by the SDSS spec- 
troscopic pipeline [^. We found 97 sources with psf g < 17.5 
or psf i < 17.5 (see Table 2) and only one is clearly RL (SDSS 
J 150324.77-1-475829.6, but it shows a spectrum with a contin- 
uum in high phase and almost missing H/3) and one is RI (SDSS 
J163323.58-I-47I858.9). One of these sources actually shows the 
broadest known FWHM H/3 (~40000 km s"\ Wang et al. 20o3) 
which exceeded the comprehension of the SDSS broad line identi- 
fier. These population A sources clearly show a small probability of 
radio-loudness. Eventual addition of such NLSyl sources (our tem- 
plate could not reduce them properly) would increase the RQ-RL 
domain occupation difference discussed in § 4. 

b) SDSS is also apparently biased against steep-spectrum 



^^ For such a task we used the SQL (Structured Query Language), instead 
of a direct selection through the Spectroscopic Query Form of SDSS; in the 
Appendix we reproduce the "where" clause we formulated. 
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Table 5. Radio Quiet and Intermediate sources in our sample that show peculiar radio morphology. 



Name 



Comments 



SDSS Jl 10717.77+080438.2 
SDSS Jl 14047.9+462204.8 

SDSS J120014.08-004638.7 
SDSS J171322.58+325627.9 
SDSS J230443.47-084108.6 

SDSS J232721.96+152437.3 



galaxy nearby; core-jet morphology 

misaligned nearby radio structure with no 

SDSS-detected counterpart 

possible core+lobes structure 

possible lobes 

RQ with elongation 

on either side of the radio core 

not in FIRST catalog; 

NVSS shows extended weak lobes 



lobe-dominated quasars JRichards et al.ll2002h . There is no reason 
to assume that any missed FRII would preferentially populate re- 
gion A in 4DE1. We tried to avoid missing RL sources with largely 
separated radio lobes and no radio detected core between them by 
carefully examining FIRST and NVSS radio maps. We consider 
that our approach is very effective in turning up all FRlIs without a 
detected radio-core at/near the position of quasar. 

c) We must also consider that SDSS does not include sources 
brighter than i ~ 15.0. This is a technical limitation imposed 
to avoid contamination of adjacent fibers by very bright sources. 
How would the omission of these bright (mostly low luminosity 
Seyferts) AGN affect our co nclusions? The latest i ncarnation of 
the 4DE1 spectroscopic Atlas JMarziani et alj|2003al) includes 215 
objects of which 61 are brighter than V=15.0 (4 RL) and 21 are 
brighter than V =14.0 (3 RL). 56/61 objects show z < 0.1 (41/61 
with z < 0.05). Most of these objects (53) show bolometric lu- 
minosity logLt,oj>44.0 erg s~^, which according to Figure 6 are 
bright enough to be RL. The RL percentage (~7%) suggests that 
SDSS exclusion of such bright AGN will not affect our conclu- 
sions. 

d) As mentioned earlier we have a total sample size of n=1770 
QSOs brighter than psf g=17.5 or psf i=17.5. We have almost com- 
plete spectroscopic coverage for all RL and RI sources in this range 
meaning that we extracted reliable 4DE1 parameters for use in Fig- 
ures 4 and 5. Similar parameters were extracted for n=333 psf g se- 
lected RQ quasars. Would inclusion of fainter psf g selected and the 
many psf i selected RQ sources change our 4DE1 definition/domain 
of RQ? We think the answer is clearly "no" for several reasons: 1) a 
random sample of 333 RQ is sufficient to define the general proper- 
ties of the RQ parent population, 2) the RQ domain defined with the 
SDSS sample is very similar to that defined from our Atlas sample 
dMarziani et al.ll2003al ; ISulentic et al.ll2007l) that shows only partial 
overlap with SDS S, 3) recent V LT spectroscopy of the H/3 region 
in high z sources JSulentic et"al . 2004, 2006; Marziani et al. 2008) 
again show the same trends as for the low redshift samples. The psf 
i selected quasars (fainter than psf g = 17.5) show too low S/N to 
allow accurate 4DE1 measures to be extracted. They are in addition 
strongly host galaxy contaminated as a class. Random examination 
of these noisy spectra give no evidence that they would change the 
general RQ properties derived from the brighter sources. 



8.3 Black Hole Mass and Radio Loudness 

The results of Figure 4 indicate that the vast majority of RL sources 
show FWHM H/3 > 4000 km s-\ with the median FWHM H/3 
of the FRII population (viewed as inclined RL sources) near 6750 
km s^^. Using FWHM H/3 and L ° measures to estimate black 

^ ^ 5100A 



hole masses and Ltoifi^Edd values implies that Ms_fr>l x 10* M© 
for RL sour ces, with a strong concentration bet ween 3x10* and 
3x10^ Mq JMarziani et ai]|2003bl ; ISulentic et al. 2006) . Edding- 
ton ratios for RL sources are restricted to hboi/LEdd < 0.15-0.30 
(see previous discussion). The RQ majority show generall y smaller 
values of Mbh and larger values of l^hoif^Edd (e.g. 



Boroson 



2002; Marziani et al. 2003b; Dunlop e t al.ll2003l ; lMcLure & Jarvis 
2004.;.Metcalf & MagliocchettL2006.) . However, some studies (e.g. 
Kording et al.l2006bn suggest that radio-loudness is not directly re- 



lated to a single variable like y\-BH or accretion rate. 

RL/RQ comparisons using spectral properties are sensitive to 
the relative contributions of Population A/B RQ sources in the 
sample under study. RQ Population B sources will show masses 
and Eddington ratios similar to RL sources while the majority of 
RQ sources (Population A) will not. The numbers quoted here 
come from analysi s of our Atlas sample (Marziani et al. 2003a|; 
ISulentic et al.l2006h . Preliminary analysis of the SDSS sample (e.g. 
Figure 4) indicate that the co nclusions wi ll be very similar to those 
summarized here (see also iLaor 1 12003| for some relevant com- 
ments). 

A full discussion an d comparison wit h other studies will 
be given in a later paper jZamfir et al.l 120081) . We emphasize the 
importance of using adequate S/N spectra and proper identifica- 
tion/subtraction of narrow H/3 for estimating black hole masses and 
Eddington ratios in AGN. 



9 CONCLUSIONS 

Three criteria have been used to isolate RL quasar samples from 
the RQ majority: 1) radio/optical flux density ratios (e.g. Rk as 
defined in § 1), 2) radio luminosity and 3) radio morphology. The 
first criterion is the least precise and was not used in this study. We 
used a combination of criteria 2 and 3 which involved determining 
the cutoff from the radio luminosities of the weakest FRII sources 
using one of the most complete RL samples ever compiled. The 
cutoff value log Li.4ghz=31.6 ergs s^^ Hz^^ agrees closely with 
the value derived in an earlier attempt using a more heterogeneous 
sample. We think this value is therefore a robust boundary for the 
classical RL phenomenon. We find many CD RL sources that are 
RL by this definition but are on the low luminosity side of most 
FRII sources. These RL CD sources cannot be FRII sources viewed 
jet-on and are either boosted RQ quasars or precursors of the RL 
phenomenon. 

We find that RQ-RL comparisons involving radio and bolo- 
metric luminosity (diagrams) yield ambiguous results about the re- 
ality of a RQ-RL dichotomy. A gap or dichotomy between the two 
populations is filled by the radio bright end of the RQ source distri- 
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bution and possible radio pre- and/or post-cursors in the zone of RI 
(and RL) sources. The optical diagnostic plane of 4DE1 provides 
much less ambiguous evidence that RL show significant structural 
and kinematic differences from the majority of RQ sources which is 
consistent with a real dichotomy. 4DE1 also shows that RI and RQ 
sources are spectroscopically indistinguishable. Our Type 1 QSO 
sample shows no evidence for an FRI population within the radio 
resolution constraints of NVSS/FIRST. 
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APPENDIX A: AN SQL SEARCH FOR QSO WITH FWHM 
H/? < 1000 KM S^ IN DR5 

We reproduce here the "where" clause we formulate to isolate ob- 
jects labeled "Galaxy" by the spectroscopic algorithm of SDSS 
(DR5); their spectra show H/3 emission line narrower than 1000 km 
s~^ in the sources rest frame. The a interval required in the query 
would translate into an observed FWHM H/3 range ~ 300-1700 km 
s~^. This takes into account the (1-l-z) scaling of line width from 
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source's frame to the observed frame, e.g. a rest frame FWHM H/3 
= 1000 km s^^ would be observed from z=0.7 as 1700 km s^^. 

SELECT . . . FROM SpecPhoto as S, 

SpecLine as L 
WHERE 

S.SpecObjID = L.SpecObjID and 

L.Lineld = 4863 and 

S.z <= 0.7 and 

L . ew > and 

L.sigma <= 11.7 and 

L.sigma >= 2.1 and 

S.psfMag_g <= 17.5 

The spectra selected this way have been visually examined and 
we kept only the bona-fide Type 1 QSO for our statistical estimates, 
as explained in the text. We selected "Galaxy" spectra based on psf 
i magnitude cut replacing the last line in the "where" clause above 
with: 

S.psfMag_i <= 17.5 



